Chaotic magnetic field reversals in turbulent dynamos 
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We present direct numerical simulations of reversals of the magnetic field generated by swirling 
flows in a spherical domain. In agreement with a recent model, we observe that coupling dipolar 
and quadrupolar magnetic modes by an asymmetric forcing of the flow generates field reversals. 
In addition, we show that this mechanism strongly depends on the value of the magnetic Prandtl 
number. 
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The generation of magnetic field by the flow of an elec- 
trically conducting fluid, i.e., the dynamo effect, has been 
mostly studied to understand the magnetic fields of plan- 
ets and stars [l[ . The Earth and the Sun provide the best 
documented examples: they both involve a spatially co- 
herent large scale component of magnetic field with well 
characterized dynamics. Earth's dipole is nearly station- 
ary on time scales much larger than the ones related to 
the flow in the liquid core, but displays random rever- 
sals. Reversals also occur for the Sun but nearly pe- 
riodically. The magnetic field changes polarity roughly 
every 11 years. Reversals have been displayed by direct 
simulations of the equations of magnetohydrodynamics 
(MHD) [3| or of mean field MHD ;3] and have been mod- 
eled using low dimensional dynamical systems 0, 0] ■ It 
has been observed recently that the magnetic field gen- 
erated by a von Karman flow of liquid sodium (VKS 
experiment) can display either periodic or random re- 
versals as well as several other dynamo regimes, all 
located in a small parameter range Q- The ability of all 
these very different dynamos to reverse polarity is their 
most striking property. This is obviously related to the 
B — > — B symmetry of the MHD equations, implying 
that if a magnetic field B is a solution, — B is another 
solution. However, this does not explain how these two 
solutions can be connected as time evolves. The VKS 
experiment has provided an interesting observation. In 
this experiment, the flow is driven in a cylindrical con- 
tainer by two counter-rotating coaxial propellers. When 
they rotate at roughly the same frequency, a magnetic 
field with a dominant dipolar component aligned with the 
axis of rotation is generated. Time dependent magnetic 
field with periodic or random reversals are observed only 
when the difference between the two rotation frequencies 
is large enough 0. We have shown that this can be re- 
lated to broken TZtt symmetry (the rotation of an angle tt 
along any axis in the mid-plane is a symmetry when the 
propellers rotate at the same frequency) Magnetic 
modes changed to their opposite by 7?.^ are dipolar ones, 
whereas quadrupolar ones are unchanged. Breaking the 



TZtt symmetry by rotating the two propellers at different 
frequencies generates additional coupling terms between 
dipolar and quadrupolar modes. If their dynamo thresh- 
old is close enough, a slightly broken symmetry is suffi- 
cient to lead to a saddle node bifurcation, which generates 
a limit cycle that connects opposite polarities. Slightly 
below this bifurcation threshold, ±B stationary solutions 
are stable but a small amount of hydrodynamic fluctua- 
tions is enough to generate random reversals. Although 
the flow in the Earth's core strongly differs from the one 
of the VKS experiment, a similar type of interaction be- 
tween two marginal dynamo modes can be considered 
and provides a simple explanation of several features of 
paleomagnetic records of Earth's reversals |13] . 

The purpose of this work is to strengthen this 
phcnomcnological scenario by displaying reversal of a 
magnetic dipole coupled with a quadrupolar mode in 
a direct numerical simulation of MHD equations. To 
wit, we consider a flow driven in a spherical geometry 
by volumic forces that mimic the motion of two co-axial 
propellers. We observe reversals of the generated 
magnetic field for a wide range of parameters when 
the propellers rotate at different speeds. We show that 
the value of the magnetic Prandtl number P,„ strongly 
affects the magnetic modes involved in the dynamics 
of reversals. Reversals that involve a coupling between 
dipole and quadrupole modes occur for Pm small enough. 
Finally, we present a minimal model for the reversal 
dynamics. 

The MHD equations are integrated in a spherical ge- 
ometry for the solenoidal velocity v and magnetic B 
fields, 

^ + (v V)v = -Vtt + j^Av + f + J-(B ■ V)B(l) 
ot flop 

^ = V X (v X B) + AB. (2) 

ot fiQa 

In the above equations, p is the density, fiQ is the mag- 
netic permeability and a is the electrical conductivity 
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FIG. 1: Magnetic field lines obtained with a symmetric forcing 
(C = 1) for i?„ = 300 and = 1. Note that the field 
involves a dipolar component with its axis aligned with the 
axis z of rotation of the propellers. 



f = /oF, where = 
, for z > 0, and opposite 



normalized 



of the fluid. The forcing is 

sin(7r s 5) , Fz — e sin(7r s c) 
for z < 0. We use polar coordinates (s, 
by the radius of the sphere a. F^ generates counter- 
rotating flows in each hemisphere, while Fz enforces a 
strong poloidal circulation. The forcing is only applied 
in the region 0.25a <| z |< 0.65a, s < sq. In the simula- 
tions presented here, sq — 0.4, b^^ = 2so and = sq. 
This forcing has been used to model the Madison experi- 
ment [ij . It is invariant by the TZtt symmetry. In order to 
break it, we consider in the present study a forcing of the 
form Cf , where C — I for z < but can be different from 
one for z > 0. This describes two propellers that counter- 
rotate at different frequencies. Although performed in a 
spherical geometry, this simulation involves a mean flow 
with a similar topology to that of the VKS experiment. 
We solve the above system of equations using the Parody 
numerical code [10]. This code was originally developped 
in the context of the geodynamo (spherical shell) and we 
have modified it to make it suitable for a full sphere. 
We use the same dimensionless numbers as in [9], the 
magnetic Reynolds number i?„i = /iocra max(|v|), and 
the magnetic Prandtl number Pm — v^qct. The kinetic 
Reynolds number is then Re = R„i/Pm- 

In a previous study with symmetric forcing (C = 1), 
we showed that different magnetic modes can be gener- 
ated depending on Pm 9]. For large enough P^, the 
dynamo onset corresponds to small Re and the flow is 
axisymmetric and generates first an equatorial dipole. In 
contrast, for P,„ small, the dynamo onset occurs when 
Re is already large and the flow involves non axisymmet- 
ric fluctuations. A magnetic field with a dominant axial 
dipole is observed (see figure [T]). In the present study, all 
the simulations are made for Re > 300, so that an axial 
dipole is obtained for symmetric forcing. 

We next break the symmetry TZt, of the forcing to 
check whether time dependent magnetic fields involv- 
ing reversals between both polarities are obtained as in 
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FIG. 2: Time recording of the axial dipolar magnetic mode 
(in black), the axial quadrupolar mode (in blue) and the equa- 
torial dipole (in red) for Rm — 300, Pm ~ 1 and C = 2.. 



the VKS experiment. Time recordings of some compo- 
nents of the magnetic field are displayed in figure [2] for 
R„i = 300, Pm = 1 and C ~ 2, which means that one 
of the propellers is spinning twice as fast as the other 
one. We observe that the axial dipolar component (in 
black) randomly reverses sign. The phases with given 
polarity are an order of magnitude longer than the dura- 
tion of a reversal that corresponds to an Ohmic diffusion 
time. The magnetic field strongly fluctuates during these 
phases because of hydrodynamic fluctuations. It also dis- 
plays excursions or aborted reversals, i.e., the dipolar 
component almost vanishes or even slightly changes sign 
but then grows again with its direction unchanged. All 
these features are observed in paleomagnetic records of 
Earth's magnetic field [ll| and also in the VKS experi- 
ment However, these simulations also display strong 
differences with the VKS experiment. The equatorial 
dipole is the mode with the largest fiuctuations whereas 
the axial quadrupolar components is an order of magni- 
tude smaller than the equatorial modes. In addition, it 
does not seem to be coupled to the axial dipolar compo- 
nent. 

Magnetic Prandtl numbers relevant to liquid metals 
are much smaller than unity (~ 10~^-10~^). While re- 
alistic values cannot be achieved owing to computational 
limitations, Pm can be decreased to values less than unity. 
We now turn to simulations using P,„ = 0.5, thus in- 
troducing a distinction between the viscous and ohmic 
timescale. The time evolution of the magnetic modes for 
Rm = 165 and C = 1.5 is represented on figure [3] (left). 
It differs significantly from the previous case {Pm — 1)- 
First of all, the quadrupole is now a significant part of 
the field, and reverses together with the axial dipole. The 
equatorial dipole remains compartively very weak and 
unessential to the dynamics. The high amount of ffuc- 
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FIG. 3: Time recordings of the axial dipole (black), the axial quadrupole (blue) and the equatorial dipole (red). 
Pm = 0.5 and C = 1.5. Right: Rm = 180, P™ = 0.5 and C = 2. 
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tuations observed in these signals points to the role of 
hydrodynamic fluctuations on the reversal mechanism. 
One could be tempted to speculate that a higher degree 
of hydrodynamic fluctuations necessarily yields a larger 
reversal rate. Such is in fact not the case. A more sen- 
sible approach could be to try to relate the amount of 
fluctuations of the magnetic modes in a phase with given 
polarity, to the frequency of reversals. Increasing Rm 
from 165 to 180 does yield larger fluctuations as shown 
in figure [3] (right) . However the reversal rate is in fact 
lowered because C was modified to C = 2. This clearly 
shows that the asymmetry parameter C plays a more im- 
portant role than the fluctuations of the magnetic field. 
For P„j = 0.5, reversals occur only in a restricted re- 
gion, 1.1 < C < 2.5, which is also a feature of the VKS 
experiment. 

Let us now investigate the detail of a polarity rever- 
sal (figure SI . Interestingly the dipolar and quadrupolar 
components do not vanish simultaneously. Instead the 
decrease of the dipole is associated with a sudden increase 
of the quadrupolar component, related with a burst of 
activity in the non axisymmetric velocity mode m — 1 
(zonal velocity) breaking the TZj^ symmetry (i.e. coupling 
the dipolar and quadrupolar families). The quadrupole 
quickly decays as the dipole recovers with a reversed po- 
larity. Immediatly after the reversal the dipole systemati- 
cally overshoots its mean value during a polarity interval. 
This behavior of the magnetic modes is typical of rever- 
sals obtained with this value of Pm and is in agreement 
with the model presented in 

These direct numerical simulations illustrate the role 
of the magnetic Prandtl number in the dynamics of re- 
versals. When Pjn is of order one, the magnetic per- 
turbations due to the advection of magnetic field lines 
by the velocity field, evolves with a time scale similar 
to the one of the velocity fluctuations. We thus expect 
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FIG. 4: Time recordings of the axial dipole (black), the axial 
quadrupole (blue) and zonal velocity (yellow) during a rever- 
sal. Rra = 165, Pm = 0.5 and C = 1.5. 



these two fields to be strongly coupled. Modification of 
the magnetic field lines due to their advection by a local 
fluctuation of the flow can then trigger a reversal of the 
fleld 14| . This type of scenario has been observed in some 
direct numerical simulations, usually performed with Pm 
of order one [l^l- When Pm is small, magnetic perturba- 
tions decay much faster and we expect only the largest 
scale magnetic modes to govern the dynamics. To illus- 
trate this argument in a more quantitative way, we have 
computed the correlation r of the most signiflcant mag- 
netic and velocity modes with the axial magnetic dipole 
for Re = 330 and 0.3 < < 1. For Pm ~ 1, all 
the modes are weakly correlated with the axial magnetic 
dipole (r < 0.3). When Pm is decreased, the correlation 
of most modes decay except the one of the zonal veloc- 
ity mode that slightly increases and the one of the axial 
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FIG. 5: Numerical integration of the amplitude equa- 
tions (|3I4I5|I . Time recording of the amplitude of the 
quadrupolar mode for fi = 0.119, v = 0.1 and T = 0.9. 



fluctuations. Thus, it is not necessary to add external 
noise to obtain random reversals. We do not claim that 
this minimal low order system fully describes the direct 
simulations presented here. For instance, in the case of 
exact counter-rotation (C = 1, i.e. F = 0), equations 
Pl4l5p do not have a stable stationary state with a dom- 
inant axial dipole. The different solutions obtained when 
H is increased cannot capture all the dynamo regimes of 
the VKS experiment or of the direct simulations when 
Rm is increased away from the threshold. Taking into 
account cubic nonlinearities provides a better descrip- 
tion of the numerical results for = 0.5. However, this 
three mode system with only quadratic nonlinearities in- 
volves the basic ingredients of the reversals observed in 
the present numerical simulations for low enough values 
of the magnetic Prandtl number. 

We thank F. Petrelis for useful discussions. Computa- 
tions were performed at CEMAG and IDRIS. 



quadrupole mode that strongly increases up to 0.8. 

We now write the simplest dynamical system that in- 
volves the three modes that display correlation in the 
low P„i simulations: the dipole I?, the quadrupole Q, 
and the zonal velocity mode V that breaks the 7^^ sym- 
metry. These modes transform as D —D, Q ^ Q and 
V under the 7^^ symmetry. Keeping nonlinear 

terms up to quadratic order, we get 



D = fiD-VQ, 
Q = -uQ + VD, 
V = T-V + QD. 



(3) 
(4) 
(5) 



A non zero value of F is related to a forcing that breaks 
the TZtt symmetry, i.e. propellers rotating at different 
speeds. The dynamical system pi4l5p with F = oc- 
curs in different hydrodynamic problems and has been 
analyzed in detail The relative signs of the coeffi- 

cients of the nonlinear terms are such that the solutions 
do not diverge when /j, > and v < 0. Their modulus 
can be taken equal to one by appropriate scalings of the 
amplitudes. The velocity mode is linearly damped and 
its coefficient can be taken equal to —1 by an appropriate 
choice of the time scale. Note that similar equations were 
obtained with a drastic truncation of the linear modes of 
MHD equations 0]. However, in that context and i' 
should be both negative and the damping of the velocity 
mode was discarded, thus modifying the dynamics. 

This system displays reversals of the magnetic modes 
D and Q for a wide range of parameters. A time record- 
ing is shown in figure [5l The mechanism for these rever- 
sals results from the interaction of the modes D and Q 
coupled by the broken 7^^ symmetry when V ^ 0. It is 
thus similar to the one described in [12I | but keeping the 
damped velocity mode into the system generates chaotic 
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